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The mass and interactions of a quantum 't Hooft-Polyakov monopole are measured nonperturba- 
tively using correlation functions in lattice Monte Carlo simulations. A method of measuring the 
form factors for interactions between the monopole and fundamental particles, such as the photon, 
is demonstrated. These quantities are potentially of experimental relevance in searches for magnetic 
monopoles. 
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I. INTRODUCTION 

Everyday experience tells us that there are no isolated 
magnetic charges, i.e., magnetic monopoles, but there 
are strong theoretical hints that they may still exist. 
Their existence would explain the quantisation of elec- 
tric charge [U, and they are an inevitable consequence 
of grand unification [1, HJ. As stable particles, magnetic 
monopoles produced in the early universe would still ex- 
ist [HQ j but there are very strong astrophysical bounds 
on their number density (as outlined by the review in 
Ref. [H). However, monopoles could also be produced 
in particle accelerators provided that their mass is low 
enough. This would clearly not be possible for grand 
unified theory monopoles, for which the mass would be 
around 10 16 GeV. It is, however, entirely consistent to 
consider monopoles which are much lighter than this, 
perhaps even around 1 TeV. They can then be produced 
in the LHC, where they are being searched for by the 
MoEDAL experiment 0. 

If the search is successful, it would open up a new win- 
dow on high energy physics. The monopoles could have 
interesting and unusual properties - such as the ability 
to catalyse baryon decay [8( - which reflect physics be- 
yond the Standard Model and yet, because they are sta- 
ble and interact strongly with the electromagnetic field, 
they would be relatively easy to study. Curiously, effec- 
tive excitations with non-zero magnetic charge exist also 
in condensed matter systems 0, [T3] , where their physics 
can be studied with simple tabletop experiments. 

To benefit from any experimental discovery of mag- 
netic monopoles and to link their properties to high en- 
ergy physics, one needs a reliable way to calculate their 
properties from theory. Calculations of their scattering 
amplitudes typically treat them as point particles and are 
hampered by their strong magnetic coupling (see Ref. (Tl| 
for a review) . However, in actual theories of high energy 
physics they usually appear as topological solitons known 
as 't Hooft-Polyakov monopoles [!,[3[, which are extended 
objects. These solutions have been studied extensively in 



classical field theory [12], but quantum results generally 
only exist in supersymmetric theories in which quantum 
corrections are straightforward. In non-supersymmetric 
theories results have been limited to leading logarith- 
mic corrections to the monopole mass at the one-loop 
level [III, [lH , although the monopole mass has also been 
calculated nonperturbativcly using numerical lattice field 
theory methods (l5| . 

In previous work [l6|, [TtJ we developed a technique to 
calculate form factors of topological solitons using lat- 
tice Monte Carlo simulations. In this paper we apply 
this technique to 't Hooft-Polyakov monopole, and calcu- 
late the form factor of the magnetic monopole for scalar 
and magnetic fields. The latter describes the interac- 
tion between the monopole and the photon, and it is 
therefore the key observable if one is ever in a position to 
study magnetic monopoles experimentally. Furthermore, 
it needs to be calculated in the full quantum theory be- 
cause the semiclassical result is that of a pointlike Dirac 
monopole, and therefore any non-trivial properties of the 
monopole appear only in quantum theory. 

We investigate the SU(2) Gcorgi-Glashow model with 
the Lagrangian 

C = —TrF^F^ + Tr[D M , $][D", $] 

-m 2 Tr$ 2 -A(Tr$ 2 ) 2 (1) 

with the covariant derivative = + igA^. The field 
$ is in the adjoint representation of the SU(2) gauge 
group and can be parameterised by the Pauli matrices as 
$ = (f) a (j a . 

In the broken phase, which occurs classically for m? < 
in the paramctcrisation chosen here, a vacuum expec- 
tation value Tr$ 2 = — m 2 /2A = v 2 forms, and the SU(2) 
symmetry is broken to U(l). In this phase, the theory has 
monopole solutions with an extended scalar field El ■ 

Even in continuum, the classical profile of this 
monopole solution must be obtained numerically except 
in the BPS limit where A — > 0. The classical mass of the 
monopole can be written as 
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where f(z) is a function of the ratio z = run/my?, mn = 
■\/2|fn| is the Higgs mass and myy = g\m\/y/X is the mass 
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FIG. 1. A two-dimensional spatial slice through the system 
with boundary conditions of Eq. ©. The magnetic field lines 
are shown, and the antiperiodicity of the magnetic field com- 
ponents can be deduced. The simulated box is shaded, with 
'image' magnetic charges shown in the neighbouring boxes. 
It can be seen that moving from the boundary on one side 
of the lattice to the other requires that the magnetic field be 
reversed. 

of the charged W^ 1 bosons. For non-zero z, the function 
f(z) needs to be calculated numerically or as a Taylor 
expansion fl8l Il9j. 

Two length scales can be associated with the quantum 
monopole. The first is the Compton wavelength deter- 
mined by the monopole mass M; the other is the core 
size of the monopole which is determined by the pertur- 
bative masses run and myv- Eq. @ shows that at weak 
coupling there is a large hierarchy between these scales, 
M > m w , m H . 



II. LATTICE IMPLEMENTATION 

In order to study the theory (fl]) using Monte Carlo 
simulations, we Wick rotate it to 4D Euclidean space 
and discretise it. Our Euclidean lattice action is 



where II + = \{1 + $) and $ = $ v /2/Tr$ 2 , giving an 
Abclian field strength tensoiQ 

2 

oi^v = -arg Tr u^(x)u v (x + (i)u^(x + v)ul(x) (5) 
and an expression for the lattice magnetic field, 

Bi = -eijkOijk- (6) 

Gauss's law for a magnetic field in standard U(l) elec- 
trodynamics is V • B = 0. In our lattice formulation, the 
corresponding equation becomes 

3 . 

[ B i ( x + *) - B i ( x )l = Pu (x) = —, (7) 

2 = 1 ^ 

where n is an integer, which can be non-zero. This means 
that the theory allows magnetic charges. In the classi- 
cal limit, these charges correspond to 't Hooft-Polyakov 
monopoles [1, Hj]. Note that even on the lattice, the mag- 
netic charge is quantised and localised in one lattice cell. 

We simulate the theory of Eq. © on a Euclidean lat- 
tice of sizc0 L 3 xT. To create nonzero magnetic charge we 
apply twisted boundary conditions on each timeslice (20j , 
while retaining periodic boundary conditions in the time 
direction. The twisted spatial boundary conditions are 

$(sc + Li) = (T l ^>{x)a l 

U(x + Li) = -a i U{x)a i ( ' 

in the ith direction, where a% is the appropriate Pauli 
matrix. These boundary conditions force the magnetic 
charge to be odd. If T is large enough then the con- 
tribution to the partition function for the simulation 
with these boundary conditions must come predomi- 
nantly from the one-charge sector as tunnelling is heavily 
suppressed. Specifically, we have a partition function 

Z tv> = 2Z (z x e- MT + ±Zfe~ 3MT + ...), (9) 



= 2^ [Tr$(x) 2 - TV$(f)[/- M (f)*(£+£)E/J(£)] 

+ \ J2 1 2 - Tr u ^(^ + m ' Tr $2 + A ( Tr $2 ) 2 - ( 3 ) 

with link matrices paramctcriscd as = I + ia a u a . Wc 
have set the lattice spacing to unity, and therefore are left 
with the gauge coupling g, bare mass m and quadratic 
coupling A as free parameters. 

In the symmetry broken phase, a residual U(l) sym- 
metry persists. We can derive link variables corre- 
sponding to this smaller gauge group [ITj, l2l| 

u„ = U + (x)U ll {x)n + (x + fi) (4) 



1 Other definitions for the effective U(l) field have been used in 
the literature [22ll . We choose this one because it preserves the 
topological nature of the original 't Hooft tensor 0] and there- 
fore defines an exactly quantised and localised magnetic charge, 
and because it is symmetric under lattice rotations. However, 
one should bear in mind that because we are dealing with an 
interacting theory, none of these expressions will be the exact 
creation operator for the real physical photon state. Therefore, 
one should really use a correlation matrix for a set of opera- 
tors with the correct quantum numbers, which is the standard 
approach in lattice mass measurements. 

2 We choose this notation to be consistent with earlier literature. 
Note that our T is not temperature, and even though we work in 
4D Euclidean space one should not interpret it as the imaginary 
time formulation of finite-temperature field theory. If one were 
to use that interpretation, the temperature of the system would 
be 1/T. 
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where 

Zi = (m,L 2 /2ttT) 3/2 (10) 

is the partition function for an isolated pointlike 
monopole (and, in fact, the usual partition function for 
a point particle at temperature 1/T - see footnote [2|) . 

Similar arguments apply to the C-periodic boundary 
conditions [23j 

$(x + Li) = <7 2 ^{x)a 2 
U(x + Li) = -<JiU{x)cfi, 

which are locally gauge equivalent to the twisted ones ((HJ , 
but not globally. These boundary conditions permit only 
even magnetic charges, with the resulting partition func- 
tion taking the form 

Z c = Z + 2Z ( y Z 2 ^e~ 2MT + ...). (12) 

Note that the magnetic field defined by Eq. (J6j> is an- 
tiperiodic with both twisted and C-periodic boundary 
conditions (see Figure [1]). In contrast, with standard pe- 
riodic boundary conditions the magnetic field is periodic, 
and therefore the total magnetic charge has to be zero. 



III. FREE ENERGY AS THE RESPONSE 
TO A TWIST 

The conventional technique for studying magnetic 
monopoles and other topological defects with lattice 
Monte Carlo simulations has consisted of measuring their 
mass via their free energy [l]| |24T - [26| . 

The mass of the monopole is obtained from the dif- 
ference in free energies in the two different topological 
sectors. This, in turn, must be obtained from the parti- 
tion functions through 

AF = F tw -F c = -ln§^. (13) 

From this, we can obtain the mass using Eqs. © and 
(fT2"| . which give 

AF = MT — In 2 — - In — + 0( e - 2MT ). (14) 

We cannot measure partition functions in Monte Carlo 
simulations (though we can, in principle, measure the 
ground state energy difference using other nonperturba- 
tive techniques [22}). Instead one can integrate along a 
path from a point in parameter space where the free en- 
ergy (and mass) of the monopole are known to vanish to 
the point where the mass is required, 



AF = 




Derivatives of the free energy along the path are given by 
expectation values, which can be measured using Monte 
Carlo simulations. Several different integration paths 
have been considered in the literature. 

We use this approach as a benchmark to compare our 
results with. We choose to vary to 2 , and integrate from 
high m 2 where the system is in the symmetric phase and 
the monopole has zero mass, to low to 2 where the system 
is in the broken phase and the monopole is massive. In 
common with most of the literature, we use finite differ- 
ences instead of a continuous derivative. The details of 
this calculation are given in Appendix [A~] 

IV. TWO-POINT FUNCTIONS 

In Ref. [l7j we introduced an alternative approach, 
which uses correlation functions calculated with twisted 
boundary conditions, and allows us to calculate not only 
the mass of the monopole but also its form factors. For 
any local operator O, one can define the corresponding 
form factor as 

/(p 2j Pi) = (P2|O(0)| Pl ), (16) 

where |p) is a quantum state with one monopole in a 
momentum eigenstate with momentum p. We normalise 
these states in a Lorentz-invariant way as 

(p'|p) = (2 7 r) 3 <5( 3 )(p'-p)F p . (17) 

The form factor is closely related to the scattering ampli- 
tude between the monopole and the particle created by 
operator O. 

In this section, we start by looking at analytical and 
semiclassical results for the form factor of the monopole 
with various quantities, then go on to generalise the re- 
sults of Ref. [17| to the present case. These results will al- 
low us to relate quantities measured in lattice simulations 
to scalar-monopolc and photon-monopolc form factors. 

A. Form factors: semiclassical results 

In the semiclassical limit, the form factor is given by 
the Fourier transform of the classical profile O c \(x) of the 
quantity O in the monopole configuration, 

/(p 2 ,Pi) = (p 2 ia(o)| Pl ) 

« M J d 3 x e l(p2 - pi) x O c i(x), (18) 

where the last line is valid in the non-relativistic limit, 
when | p»i | , | F>2 1 "C M. In this case, to which we shall 
restrict ourselves, the form factor becomes a function of 
the momentum difference k = p 2 — Pi only, as a direct 
consequence of Galilean invariance, so we will denote it 
by /(k). 
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To determine what we should expect from our lattice 
simulations, let us evaluate this for the magnetic and 
scalar field operators. First, let us take our operator to 
be 



O = Tr & 



(19) 



There is no analytic expression for the classical (non- 
BPS) 't Hooft-Polyakov monopole solution, but in con- 
tinuum the scalar field has the "hedgehog" form 



$(r) = ^=ff(r) — 
V2 r 



(20) 



where r = |x|, and H(r) is a function which approaches 
1 at r — > oo with the asymptotic behaviour [19j 



H(r) - 1 



(21) 



for m,H < 2mw- The Fourier transform of the classical 
profile Tr$ 2 = v 2 H(r) 2 has a delta function peak at 
k = 0, but otherwise it is finite and approaches a constant 
value at low momenta 



(k|Tr$ 2 |0) 



Mv 2 



H 



as 



0. 



(22) 



For more precise comparison, we will use gradient flow 
to find the classical monopole configuration $, U numer- 
ically for our chosen lattice sizes. 

Let us then consider the magnetic field B(x). This 
form factor is the most directly measurable quantity char- 
acterising magnetic monopoles, because it determines 
their scattering amplitude with photons. For a semiclas- 
sical monopole, this has the standard Coulomb form, 



^ , . 1 x 

Bx ) = -- j, 

gx A 



which has the Fourier transform 



(k|B(0)|0)=i 



AirM k 
n k 2 



(23) 



(24) 



Note that this result is the same as for a pointlikc 
monopole, which means that the semiclassical calcula- 
tion is not sensitive to the size or internal structure of 
the monopole in any way. Therefore it is not useful for 
probing the properties of magnetic monopoles, and one 
needs a quantum mechanical result instead. 



B. Form factors from two-point functions 

To calculate form factors in quantum theory, we adapt 
our method given in Ref. [l7| for obtaining the scalar 
form factor from simpler two-dimensional lattice simu- 
lations (and associated one-dimensional defects) to the 
present case. Matrix elements like / cannot be computed 



directly using Monte Carlo simulations. Instead, the ba- 
sic observable is the field correlation function, which we 
consider in the ground state |0) of the one- monopole sec- 
tor. We calculate this correlation function in momentum 
space, taking the Fourier transform in space but not in 
the Euclidean time direction, and write a spectral ex- 
pansion in terms of energy eigenstates |a) with energies 



(0(O,k)0(t,q)} 

_ v (O|6(k)|a)H0(q)|O) t(Ea _ Eo) 
4- (0|0) 



(25) 



where Eq = M is the energy of the single-monopolc 
ground state. 

Furthermore, the Euclidean spacetime is necessarily fi- 
nite in actual Monte Carlo simulations. We denote the 
length of the system in the time direction by T and in the 
three space directions by L. We apply twisted boundary 
conditions ([5]) to the spatial boundaries. In addition to 
creating an odd magnetic charge, this has the effect that 
all observables that are odd under charge conjugation 
such as B are antiperiodic, and even observables such as 
Tr $ 2 are periodic. Their momenta k and q in Eq. ([^5)1 
are therefore also quantised accordingly, 



(2ni + I)7r/L, for odd operators, 
2niir/L, for even operators, 



(26) 



with rii G Z. 

In the time direction, we impose periodic boundary 
conditions. The correlator can then be written as 



(O(0, k)0(t, q)) = |Tr U(T - t)6(q)U(t)6(k) 



where U(t) = cxp(-Ht) is the Euclidean time evolution 
operator, and Z = Tr U(T). 

With twisted boundary conditions, the states \a), \a') 
must have odd magnetic charge, and because of momen- 
tum conservation, they must also have opposite overall 
momentum k = q. The lowest such state is the single- 
particle state of a monopole with momentum k, which 
has energy 



E k = ^k 2 + M 2 ^M + —. 



(28) 



The next states in the spectrum are two-particle states 
with a monopole moving at momentum k' and a photon 
with momentum k k'. In a box of size L, the momentum 
of the photon is quantised, and therefore there is a large 
gap ~ 7r/L > k 2 /2M between the single-particle state 
and the lowest two-particle state. 
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At long time separation, we can therefore approximate Eq. (|27[) by an integral over single-particle momentum 
eigenstates |k), 

( O (0,k) O(i , q » . I / (2T % k , (2 X, (k-|<9( q) |k--)(k»|<» W |k>-^^-)-^. 

" |(2.)^'(, + k)/|^ '^- k ^;"% -^^--^-. (29) 



Similarly, we can write the partition function as 
d k i / 1 . , , . _ i , T q /* d Kj 



L 



(2^) 3 £; k . 

e v 

(2tt) 3 



k'|[/(T)|k') = E 



(2tt) 3 



2m y 



\2-7rT 



-M 



f30) 



This partition function is the individual contribution to 
the partition function from each monopolc's worldline in 
Eq. ([5]), and using Eq. (|TU)) it can be written as 



Z = Z\e 



-MT 



(31) 



To calculate the integral (f2T))) . we use the saddle point 
approximation. The saddle point ko is found by min- 
imising the action 



S(k') = E k ,(T -t) + Sk'-ki - MT 



(32) 



for given t. By approximating the integral by a Gaussian 
around the saddle point, we obtain 



/ 



(O(0,k)O(t,q)) =i(27r) 3 < 5( 3 )(q + k)x 

ri 3 fc' |/(k' - k,k')| 2 g( ko )-|( k '-k ).]vi(k )-(k'-k ) 
(2tt) 3 E u ,_ k E k , 



(33) 



where M(ko) is the Hessian matrix with components 

8 2 S(k') 



My(ko) 



dk'dK 



(34) 



k'=k„ 



In the limit of large t and T — t, the Gaussian approaches 
a delta function and we can calculate the integral 



(O(0,k)O(i,q)) = — (2^) 3 ^ 3 )(q + k) 
|/(k -k,k )| 2 1 

^ko-k^ko (27T) 3 / 2 ^(ko)' 



-S(ko) 



(35) 



3 We have corrected a typographical error in Eqs. (17) and (19) 
of Ref. [13 where the square root erroneously extends over the 
energies as well as the normalisation factor in the denominator, 
but not the action S(feo). The expressions used in the numer- 
ical analysis were correct and so the results of that paper are 
unaffected. 



where 



W(k Q ) = VdetM(k ). 



(36) 



In the non-relativistic limit k <C M, where the form 
factor is a function of the momentum difference only, we 
find 



(O(0,k)O(t,q)) 

, (2^( 3 )(k + q) (T_\ 3/2 |/(k)| 2 

L 3 \MJ E^ k E ko W(k o y 



-S(ko) 



(37) 



where we have substituted the expression (JTOj) for Z\. 

We can use Eq. (|37|) to determine the form factor from 
the field correlator. For given k and t, we obtain the 
saddle point ko by minimising Eq. (|32p . and the form 
factor is finally given by 



/(k) = ±iy/(o(o,k)o(t,-k)) 

3/4 



kE ko W(k )e 



S(k )/2 



(38) 



for O odd. The factor of i is not present for even opera- 
tors, due to parity considerations. 



C. Mass measurements 

The time separation t enters into Eq. (|3"8"|) directly 
- in parameterising the two-point function - as well as 
indirectly, in the saddle-point calculation for ko. How- 
ever, since k<Min the current calculation we can take 
the nonrelativistic limit of the action in Eq. (|32p and 
let ko = \d/T for arbitrary t. To order k, there are no 
t-dependent quantities outside of the action in our ex- 
pression for the form factor. Thus, as expected, at low 
momenta we effectively have 

(0(O,k)0(t,q)) 

= \f\ 2 c -y/M 2 +kn(t)H-y/M2 + fr-'ka(t)) 2 (T-t)+MT (gg-j 

We can use this result to conduct a fit to the corre- 
lator, noting that this is merely one contribution to the 
two-point function; the other most significant contribu- 
tion (particularly at shorter distances) will be from the 
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FIG. 2. Plot of the monopole mass obtained from the scalar 
field correlator, with twist measurement overplotted for com- 
parison. The scalar field correlator with k = (2tv/L, 0, 0) and 
permutations is used, the lowest permitted nonzero momen- 
tum. 

lightest particle that the operator O can create propa- 
gating in the bulk. Hence, if we take O = Bi then this 
will be the photon, which we will treat as massless; for 
O = Tr $ 2 it will be the bulk scalar particle, which has 
a mass mu = v2|w|. 

As always, we are assuming in using this calculation 
that the particles created by the correlation function ei- 
ther interact directly with the monopole, propagate solely 
in the bulk, or are annihilated by the vacuum. 



V. RESULTS 

Simulations were carried out using a 16 3 x 48 lattice 
with A = 0.1 and g = l/VE. For s imp licity, the pa- 
rameters are the same as those of Ref. [15| ■ For these pa- 
rameters, the theory has a second-order (or possibly very 
weakly first-order) phase transition at m 2 sa —0.27 be- 
tween the confining phase at m 2 > m 2 and the Coulomb 
phase at m 2 < m 2 . 

Configurations were created by generating a classi- 
cal 'cold' monopole and then 'heating' the configuration 
gradually towards the phase transition. Once thermalisa- 
tion of Tr $ 2 had occurred for a given parameter choice, 
the resulting configuration was used as the input for the 
next value of m 2 . In this way, a set of configurations 
was generated which could then be simulated separately 
We thermaliscd the system initially deep in the broken 
phase, and then gradually increasing the value of ro 2 , 
because moving through the phase transition in the op- 
posite direction produces extra mono pole s which would 
take a very long time to annihilate [il. |28|| . 

The system seems very susceptible to the creation 
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FIG. 3. Plot of the monopole mass obtained from the mag- 
netic field correlator, with twist measurement overplotted for 
comparison. The momenta used are k = (3n/L, tt/L, tt/L), 
and permutations. The results are similar to the scalar field 
case. 

of mctastable states, particularly long-lived monopolc- 
antimonopolc pairs as well as what appear to be excited 
states of the monopole. After thermalisation of Tr$ 2 , 
additional checks on the histogram of total charge of the 
system were carried out; fluctuations due to the finite 
volume (indicated by Eq. ) were to be expected, but 
any skewness in the distribution led us to reject the ther- 
malisation and try again. 

A. Mass measurements and comparison 

Three methods were used to measure the mass of the 
monopole. The first was the well-established response 
to a twist obtained from Eq. (fl"5j) . used previously in 
Ref. [HI and described in detail in Appendix [S] The lat- 
tices in that work were considerably smaller in the Eu- 
clidean time direction but the response to a twist mea- 
sured here is in good agreement with the L = 16 data 
that were obtained. The measurements are plotted as a 
continuous line in both Figures [5] and [3J for reasons dis- 
cussed at length in Section Mil The thickness of the line 
is the estimated error. 

It was found that, for the histograms to offer suffi- 
cient overlap that the free energy estimates (|A4[) and 
(|A5|) agreed within errors, a measurement spacing of at 
most 8m 2 = 0.001 was required. The solid line plotted 
therefore required in excess of 150 separate simulations 
to keep systematic errors at an acceptable level, although 
Eq. (|A4[) gives consistently a lower value than Eq. (|A5|) . 

The twist measurements clearly have a finite size ef- 
fect (also seen in Ref. (l5|) that affects the measure- 
ments of the monopole mass when the physical size of 
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the monopole almost fills the box. The curve of the 
twist results changes concavity as the monopole becomes 
smaller than the box size. Deeper in the broken phase 
the monopole mass behaves in a manner similar to the 
classical monopole mass. 

The errors were obtained using the methods described 
in Appendix [X] No attempt was made to account for 
the nonzero covariance between adjacent mass interval 
measurements, but it is felt that this would not give a 
major systematic contribution to the error. 

Let us now turn our attention to the use of the two- 
point correlator to calculate masses as described in Sec- 
tion |IVC] We carried out a fit to Eq. (|39)l . plus a bulk 
field which we expect will be either the scalar or the pho- 
ton, depending on the operator used: 

C(t) = Ci e -V M2 +^o(t) 2 t-^M^ + (k-k (t))^(T^t)+MT 

+ c 2 (V Bb - ik * + e - Eh ^ {T - t A . (40) 

To ensure that our error estimates are robust despite 
the clear correlations between data points at different 
separations exhibited by the two-point function, we use 
a jackknifed nonlinear least squares fit method. The error 
estimates obtained from this technique are (in the present 
work) in agreement with those from our previous use of 
bootstrapping, but there exist results demonstrating the 
robustness of the jackknifc technique for residuals that 
are not independently and identically distributed [29| . 

The length of each simulation run was about ten times 
that for each simulation used in the response to a twist 
technique discussed above. On the other hand, for the 
single point at m 2 = —0.4, the results of 150 such sim- 
ulations in two topological sectors are required for the 
twist calculation (given the conditions above of a spacing 
where the two measurements j\ and fi agree to within 
2er), whereas just one measurement in the topologically 
nontrivial sector is needed with the correlator calcula- 
tion. Added to the difficulty of thermalising every one 
of those 300 ensembles and avoiding metastability, it be- 
comes clear that it is computationally less demanding to 
use the correlator measurement deep in the broken phase 
- if it can be relied upon. Close to the phase transition, 
the finite size effects of either technique are severe in such 
a small box. We therefore leave it to future work to study 
the dynamics of quantum monopoles at strong coupling 
near the critical point [l5| . 

Based on Figure [2] it seems that there are some small 
systematic discrepancies between the twist and correla- 
tor results when the correlator of the scalar field is used. 
Given the relatively small lattice sizes used it is not in- 
conceivable that this is due to the finite size effect in one 
of the two quantities measured, but long-lived metastable 
states are another possibility. Note that we do not an- 
ticipate any major lattice artefacts playing a role in the 
monopole dynamics until ijih ~ 1, when the scalar mass 
is about the inverse lattice spacing, at which time the 
monopole will become small enough to feel the potential 
due to the discretised lattice more severely ploj ). 
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FIG. 4. Plot of the monopole form factor for the scalar field 
f<s>(k). The measurement is deep in the broken phase with 
m 2 = —0.4. For comparison, the semiclassical result is also 
shown; renormalisation conditions have been used such that 
the vacuum expectation value for the gradient flow monopole 
matches that measured in the nonperturbative simulation; its 
classical mass is then M c \ = 30.9, to be compared with M = 
34.9 ±0.1 for the quantum monopole. 

The results for the magnetic field correlator are shown 
in Figure. [3] Since the magnetic field operator couples to 
the photon, we anticipate that part of the signal in this 
case comes from a massless photon field propagating in 
the bulk. This assumption seems borne out by the failure 
of our fitting ansatz for k = (it/ L,n/ L,ir/ L), and the 
need to go to k = (3ir/L, tt/L, tt/L) to see the correlator 
expected of the monopole signal. 

Despite the apparent systematic discrepancy, the fits 
yielding the data for Figure. [3] are very good, and the 
form of the correlator given in Eq. (|3l?|) seems to be the 
right one; the long distance 'plateau' behaviour is a good 
fit. 



B. Form factor measurements 

Having studied the mass using the low-momentum cor- 
relator measurements, we now move on to the form factor 
measurements. From the results of the Monte Carlo sim- 
ulations we use Eq. (|38[) to obtain the form factor, and 
compare with semiclassical expectations. To minimise 
sources of systematic error, we use the twist results for 
the value of M in computing form factors. 

We start by looking at the scalar field form factor 
f$(k) = (k|Tr <3? 2 |0), for which there is a semiclassical 
comparison available. The classical monopole configura- 
tion can be obtained on the lattice using gradient flow. 
Following Eq. (|T51) , the form factor can then be recovered 
by Fourier transforming Tr$ 2 , for comparison with the 
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measurement from the Monte Carlo simulation. 

The minimisation was started from a classical 'hedge- 
hog' (j2TJ)) with a trivial gauge field U fl (x) = 1. A lo- 
cal minimum of the Euclidean action was obtained using 
gradient flow (see Appendix [Bl for details). The resulting 
field configurations were used to obtain the scalar field 
Tr $ 2 in the presence of the classical monopolc. As a by- 
product the classical mass was obtained (for comparison 
with the twist results above) , by looking at the difference 
in energy between topologically trivial and topologically 
nontrivial configurations, 



M c i(m 2 ) = £ t w(m 2 ) + 



4A 



(41) 



Our results deep in the broken phase are shown in 
Figure In this plot, a single value of m 2 = —0.4 
has been used for the Monte Carlo simulations, and the 
classical monopolc with the closest matching mass was 
used for the comparison. There is, unsurprisingly, good 
agreement between the two. The semiclassical agreement 
demonstrates that our technique generalises reliably from 
the relatively straightforward case of the kink to higher 
dimensions. 

The magnetic field form factor f*B(k) = (k|B|0), is 
perhaps physically more interesting. It is a vector quan- 
tity, but in continuum its direction is always parallel to 
k because of rotation invariance, and therefore only its 
length /b(&) = |fe(k)| is non-trivial. On the other hand, 
in the simulations it is easiest to consider its individual 
components [f^k)]^, but because of the boundary con- 
ditions we cannot choose the momentum to be parallel 
to a coordinate axis. Instead, we note that the length of 
the vector can be written as 



M*) = r[fe«l 



(42) 



This quantity is shown for various values of M(m 2 ) and 
k in Figure [5] For k <C mjj, we are probing wave- 
lengths longer than the monopole core size, and there- 
fore the curve approaches the expected Coulomb result 
of Eq. (J2U). In the semiclassical calculation this be- 
haviour extends to arbitrarily high momenta, which cor- 
responds to a pointlike charge, but our results show that 
in the quantum theory there is a clear deviation from the 
Coulomb result at at shorter wavelengths, when k > tuh- 
One interpretation for this is that because of quantum 
fluctuations, the magnetic charge is spread out over dis- 
tance ~ 1/ ma- 
in Figure |H1 we highlight two fixed values of m 2 and 
plot the form factor for various values of k. Changing the 
value of m 2 can be interpreted as changing the physical 
lattice spacing. Moving closer to the critical point, i.e., 
towards higher m 2 , correspond to taking the continuum 
limit. In Fig. [6j we see that closer to the continuum 
limit, the charge distribution becomes more spread out 
in physical units. On physical grounds we would expect 
that it approaches a finite continuum limit. 
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FIG. 5. Plot of the monopole form factor for the magnetic 
field fn(k), for momenta (ki, ft/L, ir/L) and permutations 
on a doubly logarithmic scale. The comparison is with the 
Coulomb result. 
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FIG. 6. Plot of the monopole form factor for the magnetic 



field fs(k), for m = —0.335 and m 
interpreted as two different lattice spacings. 
result is again included for comparison. 



-0.4, which can be 
The Coulomb 



C. A note on algorithms and performance 

Previous nonperturbative studies of topological soli- 
tons have typically employed a standard Metropolis up- 
date algorithm. It would seem, however, that excitations 
corresponding to an extended defect's worldline are not 
going to be quickly thermalised or decorrelated by up- 
dates that are local in space. Indeed, the classical topo- 
logical soliton is a solution of the field equations, so an 
obvious way to improve ergodicity would seem to be to 
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use one of the family of algorithms which relies on real- 
time dynamics. For this reason, despite the added com- 
putational cost - and the complexity arising from the 
twisted boundary conditions - it was decided to investi- 
gate the performance of a Hybrid Monte Carlo (HMC) 
algorithm. This showed promise previously when fighting 
critical slowing down in our studies of the form factors of 
critical kinks. We hypothesised that this was due to the 
fact that the defects obey the equations of motion, and so 
using an update method that integrates the equations of 
motion (or a generalisation thereof) improves ergodicity 
for observables associated with the quantum topological 
soliton. In contrast, a single Metropolis update step will 
not significantly alter the position or configuration of a 
topological defect. 

In the current situation, however, it was difficult to 
detect an advantage to using HMC. While the autocor- 
relation time was in many cases the same, the CPU time 
required to integrate a single trajectory was longer than 
a single Metropolis checkerboard sweep; the staples must 
be recalculated for every step in the trajectory. This 
poor performance may be due to our being at relatively 
weak coupling, with severe finite size effects that mask 
any critical slowing down. It may also be due to inade- 
quate tuning of the HMC algorithm to give an optimal 
acceptance rate. 

Lastly, to improve statistics for the magnetic field cor- 
relator we considered an overrelaxation step for the SU(2) 
gauge fields, coupled to an accept-reject step to account 
for the covariant derivative term in the action. We made 
use of the SU(2) move U -> UqU^Uq. This leaves the 
Wilson term unchanged when Uq = U _1 %/det V, where 
V is the 'staple' [3l|. Unfortunately we did not notice 
any substantial improvement to the statistics as a result 
of adding this step. 



probe their properties in much more detail such as scat- 
tering involving other particles. The form factor for real 
momenta, which we have calculated, is directly relevant 
for such processes. 

For the scalar we find good agreement with the semi- 
classical results, which was expected because of the weak 
coupling. The scmiclassical form factor for the photon is 
that of a pointlike magnetic charge, but our results in- 
dicate a smooth charge distribution in the full quantum 
theory. This shows that a proper quantum calculation is 
absolutely necessary in order to probe the internal struc- 
ture of monopoles using photons. The continuum limit 
deserves to be explored using the same techniques. 

It should be reiterated that although our expression for 
the magnetic field - Eq. (O - has attractive properties, it 
is not the exact creation operator for asymptotic photon 
states in the full quantum theory In principle, a numeri- 
cal approximation for the correct creation operator could 
be obtained by a diagonalisation procedure. 

In Ref . [2l[ , it was shown how to generalise the twisted 
boundary conditions to other SU(-/V)+Higgs models, N 
even. Although odd N is arguably of greater phenomeno- 
logical interest, the techniques demonstrated here should 
be equally valid in these cases. 
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VI. CONCLUSIONS 



Appendix A: Twist measurement 



We have used correlation functions to measure prop- 
erties of the 't Hooft-Polyakov monopole nonperturba- 
tively. For the monopole mass we found good agreement 
with previous studies that used the response to twisted 
boundary conditions. 

We also calculated the form factors of the monopole 
for scalars and photons. The form factor for the photon 
is physically more relevant, because it describes the in- 
teraction of the monopole with photons and the pair cre- 
ation of monopoles through the Drcll-Yan process. It is, 
therefore, the most relevant quantity for accelerator ex- 
periments It has been argued [33| that for 't Hooft- 
Polyakov monopoles it is suppressed relative to pointlike 
Dirac monopoles by many orders of magnitude. How- 
ever, in order to calculate the pair creation rate from our 
results we would have to analytically continue the form 
factor to imaginary momenta, which is not straightfor- 
ward. On the other hand, if monopoles are produced at 
the LHC, then one can envisage further experiments that 



To determine the monopole mass from the free en- 
ergy difference, as discussed in Section IIII1 one needs 
to integrate its derivative along a path in the parame- 
ter space [ljl [23 - |26j . The free energy of an ensemble is 
defined as F = — In Z, where 



VUiV<f>e~ Sl °- 



(Al) 



is the partition function. The derivative of the free energy 
difference (TLH) is therefore 



dAF 

dg 



dS ] 



lat 



dg 



as. 



lat 



dg 



(A2) 



where the subscripts indicate expectation values calcu- 
lated in the two ensembles. In our calculations, we take 
g = m 2 , yielding 



dAF 

dm 2 



[(^* 2 ) tw -(Tr* 2 > ( 



(A3) 



10 



Integrating this expression from the symmetric phase 
where M = through the phase transition to our desired 
value of m 2 will yield, in principle, yield the mass M. 
However, it is difficult to obtain reliable error estimates 
from this technique; it is important that we keep the error 
estimates under control. 

In practice, one uses finite differences instead of the 
derivative (|A3|) . The free energy difference between two 
different values of m 2 can be written in two ways, 

/ 1 = -In(e-( m 2- m ?)S.^* a ^ (A4) 

and 

/ 2 = ln( e -^- m ^* Tr * 2 ) 2 (A5) 

where the expectation values are calculated at m 2 and 
m 2 . Having established that the two measurements are 
in agreement, the change in the monopole free energy is 

AF(to 2 ) - AF{m\) = \ (jf + jf - ft - / 2 cl ) , (A6) 

where we have chosen to average f\ and f% with equal 
weights. The errors for f\ and ji for each sector are 
added in quadrature and therefore we have 

A[AF{m 2 2 ) - Af(m?)] 2 = 

7 [A/l, tw + A/| jtw + (/l,tw - /2,tw) 2 

+ A/ 2 C + Afl c + (h. c - h, c ) 2 ] (A7) 

We know M = in the symmetric phase, so we start 
summing the differences from a value of m 2 where the 
symmetry is not yet broken. Whereas the errors in each 
/ for a change from m 2 to to 2 , are themselves indepen- 
dent, there is a small nonzero covariance for two different 
adjacent steps. Some care is therefore needed when sum- 
ming all the errors in a mass measurement. 

As a check, we should make sure that our measure- 
ments of fi and / 2 are concordant, since they measure 
the same thing. This ensures that the two ensembles at 



to 2 and to 2 , are thermalised and at equilibrium, and the 
spacing is sufficiently small that the histograms of data 
for Tr $ 2 overlap adequately. Fcrrcnbcrg and Swendscn's 
work encourages us to see this process as the reweight- 
ing of a histogram of measurements, and their formula 
immediately yields 

P„ ^ 2(Tr$ 2 )e( m '" m ?) Tr<I,2 
P m 2 (Tr $ 2 ) = — 1 — s 5—— (A8) 

for the observed distribution of Tr $ 2 sampled at to 2 and 
evaluated at m 2 [Hf ; interchanging m\ and m| gives the 
expression for measurements sampled at to 2 evaluated at 
to 2 . The equality fi = f 2 then follows, but the impor- 
tance of this approach is the realisation that the mea- 
surements will not agree unless sufficient overlap of the 
histograms for Tr <E> 2 at both m 2 and to 2 , are available. 
This overlap means that the difference in the actions, 
AS = ^2 x (m 2 — m 2 )Tr$ 2 , should be relatively small. 

Less overlap means the inferred value of Tr$ 2 is an 
over (or under) estimate, being closer to the original value 
than required. This will make f\ and ji too big, so the 
free energy will be overestimated. A similar overestimate 
will occur in both topological sectors, although J2 X 
is smaller in the topologically nontrivial sector. 

Because we can use this approach to accurately in- 
terpolate the free energy of the monopole at any value 
of to 2 , the result is plotted continuously on Figures [5] 
and [3l If we encounter difficulty obtaining good statis- 
tics, then this technique would be ideal to use along- 
side parallel tempering (replica exchange Monte Carlo), 
due to the necessarily small separations between values 
of m 2 [35| ■ In our case replica exchange would have been 
usable if the measurement spacings were slightly smaller, 
but would have led to substantial wait times with our 
computer cluster. 

Appendix B: Gradient Flow 

The equations we used to minimise the classical action 
were 



4> a (x,r + St) = (f) a (x,T) + 5t 



[4(8 + to 2 ) + 8ATr$ 2 ] <p a (x,t) + \a a Uj{x,t)$>{x +j,t)U](x,t) 



(Bl) 



and 



Ui(x, t+5t) = cxp < iSta a 



Tr W a Uij(x, t)} + 2Tr{a a U. l {x, t)Q(x, t)uJ (x, t)$(x + 2, r) -h.c.j I 

staple J 



Ui(x,T). 

(B2) 



r 



The length T of the Euclidean time direction is not rele- 
vant for this process and so we could set T = 1. 
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